Topological order parameters for interacting topological insulators 



O 

(N 

o 

Q 



I 



I 

o 
o 



> 
(N 

O 
O 



Zhong Wang,''^ Xiao-Liang Qi,^'^ and Shou-Cheng Zhang^ 

'Department of Modern Physics, University of Science and Technology of China, Hefei, 230026, P. R. China 
^Department of Physics, McCullough Building, Stanford University, Stanford CA 94305-4045 
^Microsoft Research, Station Q, Elings Hall, University of California, Santa Barbara, CA 93106 

We propose a topological order parameter for interacting topological insulators, expressed in terms of the 
full Green's functions of tiie interacting system. We show that it is exactly quantized for a time reversal in- 
variant topological insulator, and it can be experimentally measured through the topological magneto-electric 
effect. This topological order parameter can be applied to both interacting and disordered systems, and used for 
determining their phase diagrams. 

PACS numbers: 73.43.-f,71.70.Ej,75.70.Tj 



Recently, topological insulators(TI) have attracted great at- 
tention in condensed matter physics|[ll-[3t]. Historically, the 
concept of the time reversal (7~) invariant TI has been de- 
veloped along two independent routes[l]. The topological 
field theory (TFT) approach was first introduced by Zhang 
and Hu[4], who constructed a microscopic model of the T- 
invariant TI in four spatial dimensions(4 + Id), and showed 
that the effective TFT is described by the 4 -H Id Chern- 
Simons(CS) actionlfsl]. In contrast to the T breaking CS ac- 
tion in 2 -H \d, the CS action in A + Id preserves the T sym- 
metry. It is now understood that this state is the fundamental 
7"-invariant TI state, from which both the three and the two 
dimensional 7"-invariant TIs can be derivedl6-8]. Using the 
simple procedure of dimensional reduction, the TFT for three 
and two dimensional 7~-invariant TIs has been constructed[6]. 
Independently, the topological band theory (TBT) was first 
developed starting from the pioneering work of Kane and 
Mele[!9], who first proposed a Z2 topological invariant within 
the non-interacting band theory. This non-interacting Z2 topo- 
logical invariant has been generalized to the three dimensional 
TIsliMl- The TFT is generally valid for systems with in- 
teractions and disorder; recently, it has been shown that the 
TFT reduces exactly to the the TBT in the non-intera cting 
limit ! 1 311 ■ The TFT has been further developed recentlv 11141- 

TIs are widely believed to be a new state of quantum mat- 
ter. Since all states of matter in Nature, including band insu- 
lators, are necessarily interacting, it is important to formulate 
a general definition of TIs which is valid for general interac- 
tions and disorder. This is especially important since an ex- 
plicit counter-example in one dimension has been constructed 
where a non-interacting topological state becomes unstable 
against interactions! 18J. After the theoretical prediction and 
experimental discoveries of the weakly interacting TIs[l-3], it 
would be most interesting to investigate topological Mott in- 
sulators where the interaction plays an essential rolei T9l - l26ll . 
In the case where strong electronic correlations play a crucial 
role, it is essential to define a general topological order param- 
eter which can determine the phase diagram of these systems. 

The TFT approach defines the generally interacting TI in 
terms of a topological term in the effective electromagnetic 
action. The angle, or equivalently the electromagnetic po- 
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FIG. 1 : Extension of G(k) to fiigher dimensional manifolds {k, u) and 
{k, u, v). (a)Extension of G(k) to G(k, u), u e [0, 1] in the definition 
of P3, with G(k, u = 0) = G(k). The 3 + Id momenta space is il- 
lustrated as horizontal Id interval. (b)Extension of G(k) to G(k, u, v) 
with G(k, u = V = 0) = G(k) in the definition of f 2. 



larization P3, can only take two discrete values in a system 
with T symmetry. In this work, we give an explicit formula 
of P3 which can be evaluated in a generally interacting sys- 
tem. Our explicit formula for the topological order parameter 
P3 satisfies the following general criteria: (1) It is well de- 
fined in the presence of interaction and disorder. (2) It takes 
quantized values, invariant against small changes of param- 
eters in the model Hamiltonian or in experiment. (3) It is 
experimentally measurable. The TBT is based on the single 
particle band states, and can not be easily generalized to in- 
teracting systems. On the other hand, the TFT is valid for 
systems with general interactions and disroder, and the quan- 
tized magneto-electric polarization can be directly measured 
experimentallvia, 1 14 l27l 12811 . The central result of this paper 
is the topological order parameter of the 3 + Id TI defined as 



Pi = ^ fdu f ^TT^^P-[Gd,G-'GdyG- 
6 Jo J (27rr 



xGdpG-^Gd^G-^GduG-^] 



(1) 



in which G = G{k, m); k - (ko, k). The momentum variables 
k - {ki,k2,k-i) are integrated over the Brillouin zone and fre- 
quency kd is integrated over (-00, -1-00). G{k, u - 0) = G{k) 
is the full imaginary time single-particle Green's function of 
the interacting system, and G{k, u) for m 9^ is a smooth ex- 
tension of G(k,u = 0) [see Fig.([Tt)], with a fixed reference 
value G{k, u - I) corresponding to Green's function of a topo- 
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logically trivial insulating state. It is convenient to choose 
G{k, u - 1) as a diagonal matrix with Gaa - (iko - A)"' for 
empty bands a and Cpp = (iko + A)"' for filled bands p, where 
A > is independent of k. Even though f 3 is a physical 
quantity in three dimensions, we introduced a Wess-Zumino- 
Witten (WZW)|29] type of extension parameter u in its defini- 
tion. Essentially, we borrowed the ideas of WZW term in co- 
ordinate space to the frequency and momentum space. Similar 
to the WZW term, we shall show that P3 is only well-defined 
modulo an integer, and it can only take quantized value of 
or 1 /2 modulo integer for an 7"-invariant insulator. The TFT 
describes the physical response of an insulator in terms of the 
eff'ective electromagnetic actionl6|] 

■^eff — '5 Maxwell + '^topo 

= j d\dt[^F,,F''^ + ^^---F,,F,,r] (2) 

where 6 can only take values of or n for a 7"-invariant insula- 
tor satisfying periodic boundary condition. We show that the 
topological order parameter P3 enters the physical response 
function by the identity: 

9 = 27rP3 (3) 

Therefore, our topological order parameter f 3 satisfies all 
three main criteria discussed above. Furthermore, the full 
Green's function entering P3 can be directly computed by 
quantum many-body techniques, such as exact numerical 
diagonalization, quantum Monte Carlo method, dynamical 
mean field theory etc. 

The fundamental T -invariant TI in 4+ Id: Even though we 
are mostly interested in the 3-1- Id physical space, the presence 
of the WZW extension parameter u hints that the topological 
structure of f 3 is inherited from the 4-1- Id TI [4]. The concept 
of a 7" invariant TI is most naturally formulated in 4 -i- Id. All 
the lower dimensional 7"-invariant TIs can be easily obtained 
from this fundamental state by a simple procedure of dimen- 
sional reduction. For this reason, we discuss the fundamental 
7"-invariant TI first in 4 4- Id. The effective TFT for the 4-Hld 
TI is given by a CS term 

5eff = ^ / d^^dt^'^^'A^dyApd^Ar (4) 

Under the time reversal transformation, A -A, Aq Aq, 
therefore, we see that this term is explicitly 7"-invariant. 

We first show that C2 is quantized to be integer for insu- 
lating system without ground state degeneracy on a 4d torus 
r"* - T12 X Suppose we penetrate a flux quanta (p - 2n 
into r34, with field strength F34 = CJ3A4 - 54A3 = (pIL^- 
On the 2d torus Tn, the 4+ Id CS term is reduced to S^^s - 
J dxidxidxae^^PA^dyAp = ^ J dxidxidxoet'^'PApdyAp 
with fj., v,p = 0, 1, 2. This is just a 2-i-ld CS term. We consider 
an adiabatic evolution of gauge potential (Ai(f), A2(f)) along 
a rectangular path C: (0, 0) (2n/L, 0) (In/L, 2n/L) 
(0, 2;r/L) — > (0,0), and we have Berry phase given by CS term 
as Ses = 27rC2. Because A/ is gauge equivalent to A/ + In/L, 



the system is actually adiabatically evolving on a torus pa- 
rameterized by (Ai,A2), Ai,A2 e [0,27r/L]. With the path C 
enclosing the entire torus surface, we have the Dirac quanti- 
zation condition C2 = integer. The only assumption in this 
argument is that we have a unique ground state; otherwise 
the Berry phase is generally non-abelian||30| and our argu- 
ment fails. From the integer quantization of C2, we reach the 
conclusion that C2 is unchanged when the the Hamiltonian is 
tuned smoothly and the energy gap remains open. 

Now we turn to the topological order parameter for such 
4 + Id TIs with general interaction. We define 

ij J {Zny 

xGdo-G-^GdrG-^] (5) 

where the partial derivative is taken with respect to the mo- 
menta kjj - (ko, ki,k2, k^, k4), and G - G(k) is the full Green's 
function. From the spectral representation of G{k), it can be 
shown that G(k) is a smooth function of k when the energy 
gap is nonzero. Under a smooth change of physical parame- 
ter, and the associated smooth change of 6G(k), the variation 
6N4 vanishes, as proven in Eq. (CI) of Ref. |6]. Therefore, N4 
is a topological invariant. Next we show that A^4 is always an 
integer. Generally, G{k) defines a map from the five dimen- 
sional k space to the space of non-singular Green's functions, 
belonging to the group GL(«, C), whose homotopy group is 
labeled by an integer: 

7r5(GL(n,C)) = Z (6) 

which is exactly A^4. Here n > 3 is the number of bands. 
Finally, we show that the identity 

C2 = N4 (7) 

which is a 4 -H Id analog of Eq.Q, holds for general in- 
teracting systems. To be specific, we consider a typical 
phase diagram shown in Fig.® for an interacting Hamilto- 
nian H - Ho(A) + H\(g), where Hq is the non-interacting 
part including terms such as tijc^cj, and Hi is the electron- 
electron(e-e) interaction part including terms such as the Hub- 
bard interaction gntintf. These two parts H(){A) and Hi(g) 
are determined by single particle parameters A - (/li, /t2, • ■ ■) 
and coupling constants g = (^1,^2. ■ " respectively. When 
(A, g) are smoothly tuned, the ground state evolves smoothly 
so long as the energy gap remains open, and therefore both 
A^4 and C2 remain unchanged, as has been discussed. Only 
when the gap closes, the full Green's function G becomes sin- 
gular, and A^4 and C2 can both change, as indicated by the 
curve ab in Fig.©. Arbitrarily picking a gapped state A in 
the phase diagram, we can find a path AC connecting A to a 
non-interacting state C, without crossing the phase boundary 
ab. Now we have A^4(A) = A^4(C) and C2(A) = C2(C). In 
addition, because C is a non-interacting state, C2 can be sim- 
ply calculated from a single Feynman diagram, which gives 
the result C2(C) = A^4(C) S |3lij33ll. Therefore, we have 
A^4(A) - C2(A) for generally interacting ground state A. 
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FIG. 2: Phase diagram on the (A,g) plane. The dark curve ab is the 
phase boundary seperating the normal insulators(NI) and topologi- 
cally nontrivial insulators(TI). The ground states are gapped every- 
where except on ab. The true parameter space is in fact infinite di- 
mensional, but this 2d diagram can illustrate the main features. This 
diagram is applicable in both 4+ id and 3-1- Id. 



Our discussion so far can be straightforwardly generalized 
to CS insulators in 2« -t- Id for integer n. In particular, we give 
here the topological order parameter N2 for the two dimen- 
sional quantum anomalous Hall insulatorl 34l 35 1: 



^2 = ? r 4^^^[^''Gd,G-'Gd,G-'GdpG-'] 
3 J (Zny 



(8) 



This formula generalizes the TKNN formulafsi*] to the inter- 
acting systems, and is equivalent to the Chern number formula 
with twisted boundary condition fsoll . It could be useful to de- 
termine the phase diagram of an quantum anomalous Hall in- 
sulator induced by interactions! 19]. This formula could also 
be used to define the chiral topological superconductor ll37l . 
where G is regarded as the Green's function of the BdG quasi- 
particles. 

The T-invariant TI in 3+ld: Now we turn attention to T- 
invariant insulators in 3H-ld. Following the logic of dimen- 
sional reduction[6], we propose a topological order parameter 
defined by Eq.([T]l. 

Now we show P3 is quantized to be integer or half-integer 
First, it can be obtained that G{k, m = 0) has the T property 
G(ko, -k, u - 0) - TG{ko, k,u- 0)"^r ' , where ^ means trans- 
position and T is the time-reversal matrix satisfying T'T = 1 
and T*T = -1. Given an extension G(k, u), with m e [0, 1], we 
consider a mirror extension, defined as G{k, u), with the prop- 
erty G(feo, k, m) = TG(ko,-k,-ufT\u e [-1, 0], for which 
Pj is defined similarly to Eq.([T]i, but with u integrated over 
[-1,0] instead of [0, 1]. It can be checked that P3 - P3 for a 
mirror pair of extensions G and G, therefore, we have 



both the same reference Green's function, so we can identify 
the manifolds (k,u = -1) and (k,u - 1) as the same manifold, 
and we have an integral over a torus in k space(fc() is still inte- 
grated over (-00, 00)). From Eq.(|9]l, we know that P3 - n/2is 
quantized to be integer or half-integer Therefore, its variation 
under a infinitesimal change of G vanishes, just like N4. It is 
well known that the WZW terms have integer ambiguity ll29ll . 
If we choose two different WZW extensions, their difference 
is generally given by the homotopy class in Eq.®, which is an 
integer. The integer ambiguity of P3 translates into the period- 
icity of under a shift of In, by the identification of Eq. (O, 
which can be obtained in similar way to Eq.©. Therefore, 
f 3 is a Z2 topological order parameter for 3 -1- Id 7"-invariant 
interacting insulators. 

Now we briefly discuss the effect of disorder. It is con- 
venient to consider a large but finite size insulator with pe- 
riodic boundary condition. To define P^,, we can use a 
twisted boundary condition|33 38] with twisted phase 0, e 
[0, 2;r], / = 1,2,3. The Green's function, as a function of 0,, 
is now a matrix whose rank is proportional to the system size 
instead of the number of bands. The P3 for disordered insula- 
tor is defined by simply replacing kj in Eq.® by 0,. It can be 
shown that this new definition of P3 reduces to Eq.© in the 
absence of disorder The analysis in the interacting insulator 
applies to the disordered insulator as well, and we have similar 
phase diagram as Fig.©, with A interpreted as the spin-orbit 
coupling and g interpreted as the disorder strength. 

Therefore, using the simple topological order parameter f 3 
expressed in term of Green's function, we have a unified pic- 
ture of TI in the presence of interaction and disorder. The 
same discussion applies to 4 -H Id and 2 + Id, where we have 
disorder-induced quantum spin Hall states, which has been 
studied in Ref.||39i|40|]. 

The T-invariant TI in 2+1 d: The topological order param- 
eter in 2 -H Id is similar to the 3 -1- Id case. The only differ- 
ence is that in 2 -1- Id we need two WZW extension param- 
eters M and V to define our topological order parameter (see 
Fig-Gh))- Given a 2 -1- Id insulating system with full Green's 
function G{k), we can extend G(k) to a 2d torus parameterized 
by (m,v); u,v e [-1,1], i.e. we define G{k,u,v) satisfying 
G(k, 0, 0) = G(k), and G(ko, k, u, v) = TG(kQ, -k, -u, -vfT'l 
The Green's functions at the boundary u - \ and v = 1 are 
fixed to be some reference value, which can be chosen to be 
trivial, say with flat bands, in the same way as in 3 -H Id. The 
topological order parameter in 2-1- Id is defined as 



2P^ = 



1 



xGdyG- ^GdpG-^Gd^G-^GduG-^] 



(9) 



where n is an integer, and Gik, u) - G{k, u) when u e [0, 1] 
while G(k,u) = G(k,u) when u 6 [-1,0]. The integral in 
Eq.® is an integer because G(k,u = 1) = G(k,u = -1) are 



P2 = 



xGdyG-^GdpG-^Gd^G-^GdrG-^ 
or 1/2 (mod integer) 



(10) 



where e^^^P""^ is the anti-symmetric tensor taking value 1 when 
the variables are ordered as (A;o, ki,k2, u, v). The cases P2 = 
and P2 - 1/2 modulo integer correspond to topologically 
trivial and nontrivial insulators in 2-i-ld, respectively. This 
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topological order parameter f 2 is valid for interacting quan- 
tum spin Hall systems in 2+ld, including states in the Mott 
regime|19]. P2 can be physically measured by the fractional 
charge at the edge of the quantum spin Hall state|41]. Analog 
of f 3 and P2 in (1 + 1)D cannot be defined because its value 
(mod integer) would vary when homotopically nonequivalent 
WZW extensions are used. 

In conclusion we have introduced topological order param- 
eters for 7~-invariant TIs in four, three and two dimensions. 
These topological order parameters are defined in terms of the 
full Green's function, and apply to both interacting and disor- 
dered systems. These order parameters take quantized values 
which are stable against small changes of physical parameters, 
and they can be measured directly in experiments. Through- 
out this paper, we assume no ground state degeneracy. New 
fractionalized topological phases can emerge when this fun- 
damental assumption is removed L42I - I44 1 . 
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